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Abstract 

& ■ 

The hyperfine effect in muon capture rate is analyzed with taking into account of the second-order 
terms in 1/M in the non-relativistic Hamiltonian. ft is shown that in the situation of the dominance of 
the squared first-order terms the interference of zero-order and second-order terms has no significant 
influence. General expression for neutrino (recoil nucleus) angular distribution in muon capture by 
nucleus with non-zero spin is proposed. High sensitivity to gp of the term related with alignment of 
the initial mesic atom is discussed in the approximation of dominating Gamow-Teller matrix element. 
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i/-) ! 1 Introduction 

O _ 

PCAC hypothesis gives a relation between the induced pseudoscalar form factor and the 
Q\ . axial vector one. In muon capture one expects gp ~ 7gA — —8.4. This prediction is in 
agreement with measurements of gp on proton and some light nuclei as 12 C and 23 Na (see, 
e.g., [l|]). However, the experiments on muon capture by 6 Li || and 28 Si |3|] yielded the 
values for gp close to zero. To clarify a situation with this suppression that partly may 
be due to medium effects further investigations are of interest. 

One of the possible ways to measure gp is the hyperfine effect (see, e.g., []!],§]). A ratio 
of muon capture rates A + /A_ from different hyperfine sublevels is rather sensitive to the 
induced pseudoscalar coupling. In particular, if Gamow-Teller matrix element dominates 
in the transition \Jn) — > \ J — In), then the ratio 

A+ _ G\ 

A_ 3(2J+1) ( r l r \ 2 ,J-l r2 ' 
j y A — 3 / + ^J~ Gp 

exhibits the high sensitivity to Gp = {E v /2M)(g P — g a — gv — 9m) or to the form factor g P . 
Here E u is the neutrino energy, and M is the nucleon mass. However, it should be noted 
that the ratio A + /A_ is proportional to (1/M) 2 , whereas the usually used description of 
the hyperfine effect is based on the non-relativistic Hamiltonian which is of first- 

order in 1/M. Thus the analysis based on the second-order Hamiltonian seems necessary 
for consistent description of the hyperfine effect. Such analysis was the first aim of this 
work. 

The second aim was to obtain the general expression for the angular distribution of 
neutrinos (recoil nuclei) in muon capture by target with non-zero spin J. Due to the 
hyperfine splitting in the initial mesic atom interference of states with different angular 
momenta F± = J ± 1/2 gives no contribution to the differential probability of muon 
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capture. Therefore, the angular distribution of neutrinos is determined by sum 



dw(n u ) x , dw F {n v ) 

where -P(-F) is the population of the state \F). Thus, irrespective of whether the process 
of capture from states with different F are detected experimentally or not, the angular 
distribution of neutrinos must be calculated separately for each state of the hyperfme 
structure. 

Explicit expressions for the asymmetry of neutrino emission along and opposite the 
direction of mesic atom polarization were obtained in Refs. 0,0,0. However, the 
ensemble of the initial mesic atoms with spin F > 1/2 may not only be polarized but 
aligned as well. This leads to the anisotropy of neutrino emission in the directions par- 
allel and orthogonal to the vector n M . The alignment effect was analyzed for transitions 
1 1/2 7r) — > 1 1/2 tt) and 1 1 7r) — > [0 7r) in a model-independent (elementary-particle) ap- 
proach in Refs.|T0, |TT| and for transitions \Jir) — > \J ± 1 7r) in the approximation of 
dominating Gamow-Teller matrix element in Ref . JT2] . It was pointed out that the align- 
ment effect is linear in 1/M and is very sensitive to gp. Nevertheless, the contributions 
of non-leading matrix elements, in particular, of velocity terms, are evidently of great 
interest. 



2 Second-order Hamiltonian 

In the initial ls-state of mesic atom a captured muon is described by the 4-component 
wave function. It can be written in the non-relativistic approximation as (h = c = 1) 

MW»t) = V*.(r M ) ( ) e***, (3) 

where is the total muon energy, including its binding energy in the atom, ip^a^) is 
the two-component spinor, and a M is the projection of the spin s M = 1/2 onto the z axis. 
A final neutrino with momentum \l v and projection <j v of the spin s u = 1/2 onto the z 
axis is described by the wave function 

VvK, r„, t) = u„(k v , a u )e^- EU \ (4) 

where u u (k u ,a v ) is a 4-spinor, and E v = k u is the neutrino energy. Assuming that the 
weak nucleon-lepton interaction is pointlike, we introduce the lepton current acting on 
the nucleon 

Jx(a„cr u )e- l{kM ^ = r, t) 74 (l + 7s )^K, r, t), (5) 

where u = E^ — E v . Thus the relativistic Hamiltonian for the nucleon into a nucleus 
involving in capture of muon being in the spin state |cr M ) with emission of neutrino in the 
state |k„, o v ) is 

H = M/3 + op + U(r) + H w , (6) 
Hw = ^^jx(*», a v )e-^ + ^iir x )U (7) 

where G is the weak- interaction coupling constant, 8q< is the Cabibbo angle. The lowering 
operator f_ acts in isospin space and transforms a proton into a neutron. For simplicity 
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we take the nucleon-nucleus potential U(r) in a central form. The operator of the weak 
nucleon current is given by 

r ( . 9m i -9Pj 

Ta = 74 ( 5V7A + ^[ a ^P k P ~ 9AJXJ5 ~ l ~ k ^ 

It involves the matrices o\ p = (7a7 p — / ~f p '~fx)/2i, the muon mass m, and the 4-momentum 
transfer 

k\ = i>\ - fix = (k„, -iu), (9) 

where z/a and fix are the 4-momenta of the muon and neutrino, respectively. The form 
factors of vector interaction gy, axial- vector interaction g A , weak magnetism g^, and 
induced pseudoscalar interaction g P depend on k 2 = k\k\. 

To go to the non-relativistic description of intranuclear nucleon we perform the Foldy- 
Wouthuysen transformation of relativistic Hamiltonian. Keeping the terms up to second- 
order in 1/M we obtain for non-relativistic Hamiltonian for j-th nucleon the following 
expression 

if 1 1 : , 



H = M+ 3 



2M 



U{r ) + 



8M 2 



AM 



G cos 6, 



V2 



C_ -ik„r,- 



G' v + G' P {trjTL v ) + 



+ iG l {n v [(T j x ^]) + G 2 (n^) + ^p^ApA (^n,) 



+ 



G >j + ^ff + l 9il<7j x ff ] + 



#2 



[cr X M ] X M 



C/'fo) 
~A~M 2 ~ 



where n u = k u /k u , iij = Yj/rj, and 



(10) 
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G' v = g v 1 + e - 



AM 



£\- 9m 



2M 



G' A = g A ll- £ -\-e [g v {l - |) + ^(1 - 77) 



G' P = e ( fop - </ A - #y)(l - |) - g M (l -v)\, 



9y — 9v - 29a, 



g A = g A (l + -) + g P -, 



(12) 



Gi = ^(gv + 9a + 29m), 



G2 — ~2~9a, 



9i = 2^ £ 9a ~ v(9v + 2g M )), 



2M' 
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2M' 



9 a, 
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-,9a, 



(13) 



The non-relativistic Hamiltonian of muon capture by free nucleon (i.e. U (r) = 0) with 
second-order corrections was earlier obtained in Ref. 1T3H0 as well as in the other forms in 



Refs. 14, 11 



Our result differs by the factor 1/4 in the term m/(4M)e in the definition of G' v and by the sign of 33. 
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The zero-order and first-order terms are proportional to the operators 1, <Tj, pj, and 
(a-jpj). The second-order terms involve the additional set of operators: [crj x pj], [[<Tj x 

Pi] x Pj], CT i( n ^Pi) 5 (°i n j), and [°j x 



3 Neutrino angular distribution 

Let \JfMf) be the wave function that describes the internal state of the final nucleus with 
spin Jf and its projection Mf onto the z axis. At the same time, the initial state of the 
mesic atom with total angular momentum F involving a nucleus with spin J, and a muon 
is represented as 

\F) = E C&^JJMWM- ( 14 ) 

£ Mia,, 

Polarization and alignment of the ensemble of the mesic atoms with given total angular 
momentum F are described by density matrix 



p^(F) = a e (F)aUF), 



or by spin-tensors 

r Qq (F) = £ C^ QqPce (F), r 00 (F) = 1. 
«' 

The energy of the neutrino emitted in the fixed transition \J t ) — > | Jj) is 

1/2 



(15) 



(16) 



£7„ = M 



1 + 



2Q, 



M 



1 



2M 



+ . . . 



(17) 



where Q M = Mj + — My, Mj and My are the masses of the initial and final nuclei, 
respectively. Using the Fermi rule, we find that the differential probability of muon capture 
per unit time from a given state \F) of the hyperfine structure has the form 

dw F (n,,) 1 



x E 

a v M s 



(2nf 



El 



l + E u /M } - 



It is convenient to introduce the multipole expansions for the matrix elements of all 
operators appearing in the non-relativistic Hamiltonian. Generalizing the definitions of 
Ref. II for the operators of the first-order Hamiltonian we get for scalar operators 



(4vr) 3 / 2 E(-l)"^ m (n,)^ 



i^jPj) 
iU'ir^ia-jiij) 



f 

um 



[Quu] 
[Ouu,p] 
[Ouu, r] 



(19) 
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and for g-th spherical components of vector operators 



(J f M f \J2e- 
i=i 



a 31 

Vn 

A°3 X Pi] q 
[[tTj X Pi ] X pj] q 

U'trJltTj x nj} q 



(4vr) 3 / 2 



'12 nw* (n \ srr uM r 1 



aM n J f M ! 

uM 



(20) 



[lum] 

[l-ura, ap] 
[lwu, ap 2 } 

[lwu, r] 
{lwu, up} t 

The quantity {lwu, up} is linear combination of the reduced matrix elements [kwu, ap) 
(k = 0, 1 or 2) 



{iwu,a P } = (-i)«-« EC-^Ve 



2fc + 1 



2m + 1 



x 



X (-/w [/(ui W-l lfc, ltt)[fc W-l M, fjp] 



(21) 



— + 1 {7(w ty+l lu)[kw + l u, ap 



defined by 



CjiM!uM[kwu,ap] = V2{J 3 M_ 



3=1 



f\\ 4vr ^Jr»\.™*3 



X 



A? 



:22i 



Note that [Ouu, ap] = —y/2[0uu,p\. We use the normalized Racah function U(abcd, ef) 



(2e + 1)(2/ + 1) W(abcd, ef). All reduced matrix elements are real- valued quantities, 
provided that the nuclear wave functions are transformed under time reversal in standard 
way |T5| 

f\JM) = (-1) J+M | J - M). (23) 

Let be the unit vector along the axis z of orientation of the ensemble of the initial 
mesic atoms with given spin F. Polarization and alignment of this ensemble are fixed by 
spin-tensors t w (F) and t 2 q(F) or normalized parameters 



Pi(F) 



lF + 1 



r 10 (F), p 2 {F) 



(F+l)(2F + 3) 
\ F(2F-1) 



T20{F). 



(24) 



Both parameters are equal unity if only the state corresponding to the maximal projection 
£ = F is populated (pi(F) = (£,)/F is the conventional polarization). 

After straightforward calculation we get from Eq.flTBD the general expression for neu- 
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trino angular distribution 



x < 



U{hiFK, FJ,) Y^U(J f u'J z K, J t u) C$> K0 V(u)V(u') - 



^C^ K0 EC: %A(wu)V(u') + 



+ U-lf^Y^^Y'^CwQKoUilw'uK, u'w)A{wu)A{w'u') + 



+ ^ ]T JH^N +l)(2u' + l)(2w + l) < 



ww'N 



w' v! 1 
w u 1 
JV X 1 



k /WO /-<iV0 v 

' °wOAfo°ioft:o A 



x A(wm)A(wV) I + 

' Ji F 1/2 ' 

+ E\/ 6 ( 2iV + 1 )( 2J * + X )( 2F + X ) 1 Ji F X / 2 I S U(Jfll'JiN, JiU) x 

N K 1 uu ' 



N 



C%° 0l0 C:»° m V(u)V(u>) - -^J2C^ K0 U(ulu>K,wN)A(wu)V(u') 
~ 2,[lYA(wu)V(u') Y,C^ C$? 0W U(ulu'A,wN)U(llNK, 1A) - 

V w A 

- ^(-l)"'-"^ %^(-l)»'-Ti [/(ltt;W,^)A(H^(«'V) + 



/9 . I u w 1 

+ Y E^oW 3 ^ + l)(2u' + l)(2w + 1) «' 1 
■ # if 1 



> x 



x A(wu)A(w'u') 



> . 



(25) 



It is a series in the Legendre polynomials Pq(cos9) = 1, Pi(cos#) = cos8, P2(cos8) 
3(cos 2 9 — l)/2. . . , where 9 is the angle between n u and n M . The constant is 



C^ = 8(Gcos9 c y 



mZe 



2 \ 3 



R{Z)El 



(26) 



l + m/Mj l + E v /M/ 

where Z is the charge of the initial nucleus, and R(Z) is the correction factor for its 
non-pointlikeness. The amplitudes V(u) and A(wu) for the second-order Hamiltonian are 
given by 

G' V [0UU] + G 1 {1U :° P} + G 2 ^fi, if TT^-l)" = 7T f , 



V(u) 



G' P {lu} + g' A 



, \0uu,p] [Omm, r 



M 



(27) 
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m u i i [lwu, ap 2 ] . {lwu,ap} ., , 1 , w 
G' A [lwu]+g 2 l m' 1 + 9 * — jvf 1 ^ if vr 4 (-l) w = 7T/, 

A(wu) ^ ^ [lwu,p] [lwu,ap] [lwu,r] l _S ' 



where 



9v M +91 — Id + 9v AM 2 ' if n ^~^ W = _7r /' 



{1m} 1 / f [1m— 1m 



{1m, P } > = ~\ 3 ( 2u+ l) \ l lu - lu >P] } + 
{lu,cxp} j ' I [lu— lu, ap] 



, /( , ( [1«+1«] 

V^« + ij I [i u+ i U)(Tp ] 



(29) 



7Tj and 7Tj are parities of initial and final nuclear states, respectively. 

The case of transition Jj — > Jf = Jj±l, 7Tj = 7T/ with dominating Gamow- Teller matrix 
element [101] was considered in Ref . [ 0| . The neutrino angular distribution has the form 



dw F {n v ) _ C At [101] 2 2J / + l fl B f ^ fn n s , 
dn ~ 12vr 2 X + 1 (a ° + aiM^X*^) + 



+ a 2 (F)p 2 (F) 



3(n„iV 2 



(30) 



where 



a = G'^ - \g' p )C x (J u J f , F), (31) 

ai = —G A C 2 (Ji, Jf, F) — G' A (G' A — G' P )Cz{Ji, Jf-, F), (32) 

0-2 — G A G P Ci(Ji, Jf, F). (33) 

The explicit expressions for coefficients Gj(Jj, Jf, F) are presented in Ref. ||12||. Due to 
high sensitivity of the anisotropy a 2 to G' P the form factor g P can be measured by this 
way. Note that this effect is of first-order type {G' A G' P ~ s). Using the general expression 
(p5|) it is easy to estimate the contribution of non-leading matrix elements for any given 
transition. 

4 Hyperfine effect 

The rate of muon capture from hyperfine state \F±) (F± = Ji ± 1/2) of mesic atom is 
given by isotropic term of differential probability ( p5|) , i.e. 

It can be represented in the form proposed in Ref. || 

A F = A + 5A F , (35) 
where the statistically-averaged muon capture rate is given by 

x= ^§tt^( x2(m)+z/2(m) )' (36) 
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and the hyperfine increment takes the form 

SA - C » 2J f ± 1 J ^ + X ) + 3 / 4 ~ F ( F + X ) 



x 



X 



•' 2 2J { + 1 2J i {J i + 1) 

E ((•*(•* + 1) + «(« + 1) - J/( J/ + 1)) - ^4) - 

- 2^(J,+ J f +u + 2)(Ji- J f +u + l){Ji + J f -u)(J f - Ji+u + l) x 



m + 1 



where 



x(u) 



2u 



2u+ 1 



V(u) 



2(u + l) 
\ 3(2«+l) 



A(mm) + \ ~A(u - 1m) 



I2(u+ 1) 



^(m 



2u 



A, 



2 2Ji + 1 2(J, + 1) 



3(2m + 1) 
■ be brc 

+ 1) ^M+l 



A (mm) 



~A(m + 1m) 



2m + 1 

However, the muon capture rates A± may be brought also to the form 
C M 2J f + 1 



x 



(Ji + J f -u)(J f - Ji + u + l)x(u + 1) + 



+ J(Ji + Jf + u + 2)(Ji-J f + u+ l)y(u 



_C ll 2J f + \ 1 v 1 
2 2J, + 1 2Ji^-M + l 



x 



Putting 



x y (Ji + J f + u + 2)(Ji -J f + u+ l)x(u + 1) 



- \J(Ji + J f -u)(J f -J i + u+ i)y(u))' 
\ M u (u), if n i (-l) u = 7r f , 

[ ~M U (-U), if 7ti(-l) U = -7T/, 

M u {-U-1), if 7T 4 (-1)" = 7T/, 
M U (U+1), if TTsC-l)" = -7T/. 



X M 



(37) 



(38) 
(39) 



(40) 



(41) 



(42) 



(43) 



we get the formulas || for muon capture rates from different hyperfine states expressed 
in terms of the amplitudes M u (k). 

On the other hand, using Eqs.(|2"TD, ( [ZED we obtain by this way the explicit expressions 
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for the amplitudes M u (k) with the accuracy up to the second-order terms in 1/M 



+ 



I2u + 1 



, u ^\ \\u-lu,p\ , / u(u + 1) n [lu + lu,p] 
9v - 2^TI G V M + V 3(2n + 1) G2 M 

'2u + 1 / u n u+1 \ [lu — 1 u, up] 

[2^TT Gl ~ 9l ~ 2(2u + l)^ 3 ) M + 

/ u(u+l) ( 1 \ [lu+lu,<rp1 
+ V 3(2^ + 1) r 1 + 2^ 3 J M + 



1 / (u- l)(u + l) [2u-lu,qp] _ 1 / (u + 1)(u + 2) [2u + lu,qp] 
+ 2V 3(2u + 1) 93 M 2 V 3(2u + l) ^ 3 M 



u + 1 \luu,ap 2 } 2u + 1 [lu — lu, 
"^~^ 2 M 2 + V^ _ ^ 4M 2 



M u (-u - 1) = ^2^1 |v^+TG" y [0H + ^G' A [luu] - 

\2u + 1 / , -u + 1 ^ \ [1-u + l u,p] / u(u + 1) n [lu—lu,p] 

+ l2un Gl " 2(2u + l)^ 3 J M 

/u(u+l) ( n . 1 \ [lu—lu, crp] 
" V 3(kTlT l Gl + 2^ 3 J 14 ~ 

1 \ u{u-T) [2u-lu,crp] 1 / u(u + 2) [2u + lu,ap] 
2\/3(2u+l)^ 3 M + 2V/3(2u+l)^ 3 M + 

lu \luu,ap 2 } I2u + 1 [lu+lu, rll 
+ V3 * JP - V~3-^ 4Af a J • 

/ (u- l)(2u + l) ~ [2u-2u,ap] 1 \ {u + l)(2u + 3) [guu, <rp] 
+ Y 6(2u - 1) 93 M 6y 2u-l 93 M 

2u+l [lu — lu,ap 2 } i— [Ouu, r] u + 1 [luu, rll 
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M > + « = VSiTTl V^T^ - 5Tl^) l lM+1 "l 




ff, 4 + \/mTT (si - J*) ^ + 



u , [1mm, p] [u ( 1 \ [lw^pj 



3 



^ M + \3{ 9l + 2 93 M ^ n 



(m + 2)(2m + 1) [2m+2m,ctp1 1 u(2u - 1) [2mm, <rp] 
6(2m + 3) ^ 3 M + 6V 2tt + 3 93 M 



2u+l \lu+lu,ap 2 ] i — pr [Omm, r] /m [1mm, rl | 
"V^^ 2 M^ + V ^ TT ^ 4M^ + \/3^ -JM^Y 

We see that the second-order terms do not change the general expressions (35)- (37) or 



ETD for muon capture rates from hyperfine sublevels but only are added to the 
known zero-order and first-order terms in the amplitudes M u (k). 

To study the magnitude of interference between the zero-order and second-order terms 
with respect to the squared first-order terms let us consider as an example a transition 
3/2 — ► 1/2 between the states with the same parity The muon capture rates are given by 

A + = ^ (m 2 (2) + -|=M 1 (2)M 2 (2) + 1m|(2) + ^Mf(-3)) , (48) 



A_ = HSfi (m 2 (-1) + 1m 2 (2) - ^Mx(2)M a (2) + ^M 2 (2)j , (49) 

The first-order term ~ G" P [101] may exceed in the amplitude Mi (2) both the zero-order 
term ~ (G' A — 2G' P /3)[121] and the other first-order velocity terms due to the dominance 
of Gamow-Teller matrix element [101], on the one hand, and the high absolute value of 
the form factor gp, on the other hand. Just in this case the squared first-order term 
~ Gp appears in the numerator of the ratio (H). Then it is apparent that the significant 
influence on the magnitude of the squared amplitude Mi (2) is beyond the capabilities of 
any second-order term in it, as 

M x (2) = M-^ G' P [101] - (G' A - lCp) [121] + 



+ V2^g A 3 g 3 J M + ^9v M "+ 



(50) 

1 ( a | 1 yill,gp ■ 1 _ 211, gp 3 [231, op] 
[121, ap 2 } m [011,r] 1 [lll,r]l 

It is clear that the contribution of second-order terms does not exceed several percents. 
However, the terms associated with nucleon-nucleus potential, which were never consid- 
ered earlier, deserve special attention. 
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5 Summary 



The second-order in 1/M non-relativistic Hamiltonian of muon capture is obtained. For 
the first time the terms associated with nucleon-nucleus potential are considered. General 
treatment of the hyperfine effect in muon capture rate is given with taking into account 
of the second-order terms. In particular, the generalized expressions for the usually used 
amplitudes M u {k) are presented. It is shown that in the situation of the dominance of 
the squared first-order terms the interference of zero-order and second-order terms has 
no significant influence. However, the second-order corrections may be of interest due to 
their sensitivity to the nuclear potential. 

General expression for neutrino (recoil nucleus) angular distribution in muon capture 
by nucleus with non-zero spin is proposed. In the approximation of dominating Gamow- 
Teller matrix element we discuss the high sensitivity to gp of the term related with 
alignment of the initial mesic atom. Such alignment may be realized in capture of muons 
by oriented atoms. Using the general formula for the angular distribution the contribution 
of non-leading matrix elements to alignment effect can be easy estimated for any given 
transition. 

I am grateful to Prof. J.P.Deutsch for his interest to this work and useful discus- 
sions. The work was supported in part by grants 96-15-96548 and 96-02-17517 of Russian 
Foundation for Basic Research. 
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